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Dipolar interaction between the magnetic moments of electrons is studied as a source for electron 
spin decay in quantum dots or arrays of quantum dots. This magnetic interaction will govern spin 
decay, after other sources, such as the coupling to nuclear spins or spin orbit coupling, have been 
eliminated by a suitable sample design. Electron-electron (Coulomb) interactions, important for 
magnetic properties, are included. Decomposing the dipolar operator according to the symmetric 
group of electron permutations allows one to deduce vanishing decay channels as a function of 
electron number and spatial symmetries of the quantum dot(s). Moreover, we incorporate the 
possibility of rapid phonon induced spin conserving transitions which crucially affect the temperature 
dependence of spin decay rates. An interesting result is that a sharp increase of the spin decay rate 
occurs already at relatively low temperatures. 
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I. INTRODUCTION 



One of the proposals to realize qubits for quantum 
computingi uses the electron's spin in semiconducting 
solid state nanostructuresiSi^ This approach could benefit 
from traditional electronics device experience and allows 
for straightforward scalability. Recent efforts have suc- 
ceeded to demonstrate controlled preparation and detec- 
tion of singl^ and of pairs^ of electron spins in quantum 
dots. Aiming for long running coherent computations, 
however, solid state based devices somewhat suffer from 
relatively short spin decoherence and relaxation times as 
a drawback, for example, compared to nuclear spin based 
qubits.^ Exponential decay of the upper Zeeman level 
population has been observed in GaAs quantum dots'* 
over times not exceeding 10~^ sec. Therefore, it is im- 
portant to know (and, if possible, to control) any kind 
of mechanism causing spin relaxation in solids, and par- 
ticularly in semiconductors. Previous theoretical work 
has valued several contributions. They can be subdi- 
vided into two classes regarding the magnetic mechanism 
to mix spin states as a source for spin decay: (i) spin- 
orbit coupling piSiSiiS also, recently, in its interplay with 
the electron-electron interactiou j^^i^^ or (ii) coupling to 
nuclear spins. The latter can act through the spin-flip 
Overhauser effect by hyperfine interactionsii24i^*i^ii& In 
principle, both magnetic sources for spin mixing can be 
eliminated by a proper device design and by the choice 
of the semiconductor material. Non-vanishing spin-orbit 
splitting can have several causes in semiconductors which 
we briefly address: p-type bands, i.e. usually valence 
bands, may split by spin-orbit effects, arising near the 
nuclei. Secondly does the lack of spatial inversion sym- 
metry produce spin splitting even of s-type bands, ci- 
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ther by the Dresselhaus mechanism arising in the absence 
of crystallographic centrosymmetry as in Zincblende or 
Wurtzite structures. The latter particularly refers to all 
III-V semiconductors with GaAs being the most striking 
example. Also devices lacking structural inversion sym- 
metry, e.g. near surfaces or in asymmetric quantum wells, 
producing internal electric fields show spin splitting due 
to the Rashba mechanism. Fortunately, this latter spin- 
orbit source may be suppressed by fine tuning suitable 
gate voltagesiii The goal to avoid spin-orbit coupling ef- 
fects therefore suggests using spins of conduction band 
electrons in Si or in Ge in carefully symmetrically pre- 
pared structures. Also the attempt to avoid coupling to 
nuclear spins favors the use of Si or Ge: their natural iso- 
topic mixture contains nuclear spin / = to more than 
95% (Si) or more than 92% (Ge), respectively^^ 

With this work we consider the effect of dipolar inter- 
actions between the magnetic moments of electron spins 
which, for fundamental reasons, cannot be removed by 
design. While considerably weakerJA than the above 
quoted mechanisms this interaction unavoidably causes 
spin relaxation and, in the absence of other magnetic in- 
teractions, combined with the never vanishing electron- 
phonon coupling(2S will set the ultimate limit for long 
time quantum computations using electron spins, even 
in optimally designed structures. We study transitions 
between energy levels differing in their total spins^-'^ and 
disregard here effects associated with transitions between 
Zeeman levels (which conserve the symmetry of many 
body electron levels, see below) at finite magnetiza- 
tions when an external magnetic field is applied. Af- 
ter introducing the model in Sect. ^1 we reveal circum- 
stances of particular spin stability w.r.t. dipolar interac- 
tions (Sect. lTTL^ . also regarding excited (many electron) 
states in Sect. 



nil Dl depending on the electron number 
and on the symmetry of the single or the ensemble of 
quantum dots. 

We explicitly include Coulomb interactionsiiiiS due to 
their importance for magnetic properties. For example, 
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they can cause total ground state spins greate r^^i^^i^^^i^^ 
than S" = or 5 = 1/2, as expected for even or odd 
numbers N of non-interacting electrons. Here, we take 
the electrons to be confined inside one quantum dot or in 
different quantum dots. As a complementary approach to 
the non-interacting or weakly interacting regime we focus 
on strong Coulomb interactions at low electron densities, 
where pocket stateS)2SiSL2& offer a reliable description of 
many body states localized by Coulomb repulsion, even 
in single quantum dotsi^ Pocket states are briefly re- 
viewed in Sect. lIIl'ni for the present purpose to determine 
matrix elements of the dipolar interaction in Sect. IIII'DI 
Pursuing most of the foregoing theoretical work on 
electron spin decay in quantum dots we consider in 
Sect. IIV Al phonons (which themselves cannot change 
spin states) to provide the transition energy between dis- 
crete dot levels. Contrary to extended bulk situations, 
this transition energy is much bigger than mere mag- 
netic energies which is is one reason for the relative sta- 
bility of quantum dot compared to bulk electron states, 
in accordance with experimental facti^^ Spin changing 
transitions due to the combined action of dipolar energy 
and phonons are discussed in Sect. IIV Bl Generalizing 
previous results we account for rapid spin conserving ex- 
citations of the electron system induced by phonons that 
occur already at relatively low temperatures; these tran- 
sitions turn out to govern predominantly the temperature 
dependence of spin decay times, discussed in Sect. IIV CI 
Finally, we summarize and value our findings in Sect.lVl 



II. MODEL 

Specifically, we consider the iV-electron system 

N 



i=l ^ ^ i^j 

confined by the potential v{r) which is supposed to de- 
scribe a single quantum dot or more complex situations 
of many quantum dots, such as for example N quantum 
dots each containing a single electron. To be realistic, 
particularly regarding magnetic properties, we include 
interactions between electrons w{r) = 7-?""^, depending 
on the static dielectric constant k of the host material; 
Coulomb interactions are always considerably stronger 
than dipolar energies. Moderate screening, not reducing 
the interaction range to values smaller than the electron 
separation, will not affect qualitatively our results. In 
Eq. J^l, Pi and are d-component momentum and po- 
sition vectors, depending on the dimensionality d of the 
quantum dot wave functions (in heterostructurcs, d ^ 2); 
m* is the band electron mass. 

Notice that at strong Coulomb repulsion, which is the 
focus of this work, N electrons Wigner crystallized22i2£ 
in a single quantum dot become in their theoretical treat- 
ment at low energies very similar to the case of TV elec- 
trons localized in separate quantum dots. The essential 



physics of both situations is captured by an antiferromag- 
netic Heisenberg lattice model,^§^ Eigenstates of Hq ex- 
hibit well defined spins S and can be classified according 
to the eigenvalues of the z-component Sz and the square 
of the total spin operator S — X]i!=i yielding eigen- 
values Sz and S{S + 1), respectively. With SU(2) sym- 
metry in spin space, Zeeman multiplets —S<Sz^ +S 
are degenerate. We index eigenstates |(/)„> and eigenval- 
ues En of Hq by n, taken to incorporate the values of 
S and Sz- Transitions between |0„> and \4>n'> may or 
may not change S. The present work focuses on inelas- 
tic transitions that change the total spin values S S' 
rather than on transitions within a Zeeman multiplet. As 
already mentioned, many-electron ground states \4>n=o> 
may exhibit total spin values So > 1/2 as a result of 
electron-electron interactionsi22i2^iSii2^ 



III. TRANSITION MATRIX ELEMENTS 

In order to satisfy the Pauli principle an iV-fermion 
state <ri, si, . . . , r^r, sat | must belong to the A2 = 
[1^] representation of the symmetric (permutational) 
group Sn with respect to permutations p G Sn of the 
particle enumeration, {1,...,A^} {p{l) , . . . , p{N)} , 
see Ref. l32l When permuting only spin coordinates 
{si, . . . , Sn} {Sp(i), . . • , Sp(Ar)} the state |0„> trans- 
forms according to the irreducible representation (par- 
tition) F = [iV/2 + S,N/2 - S] of Sn for spin- 
i Fermions2LSa at given S = | 1 , . . . , 7V/2 for 

{ ckM } Correspondingly, when permuting only posi- 
tions {ri, . . . ,rAr} -> {rp(i), . . . ,rp(jv)}, |0«> transforms 
according to f = [2^/2--S^ issj (^j^i^ F x f containing the 
A2-representation). We notice that total spin changing 
transitions require altering the wave function's symme- 
try, which necessitates operators acting simultaneously 
in position and in spin space (by contrast, transitions 
within a Zeeman multiplet leave unaltered the symme- 
tries of wave functions). 



A. Dipolar energy 

Here, we investigate the dipolar interaction be- 
tween electrons. As seen in Eq. Q below, it contains 
products of position and spin operators and, indeed, 
mixes spin states. However, it is by far too weak to pro- 
vide the energy 

separating quantum dot eigenlevels. Focusing on Si, 
we consider in Section llVI acoustical deformation poten- 
tial phonons'''^ to supply the necessary transition energy. 
Unaided electron-phonon coupling, though, does not mix 
spin states and thus leaves spins unaltered. Eventually, 
it turns out that dipolar interaction, as a result of its 
smallness, ensues considerably smaller transition rates at 
low temperatures than, for instance, nuclear spin induced 
spin mixingji^ 
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The operator of the dipolar energy 



(2) 



is, as required for identical particles, invariant with re- 
spect to permuting the electron enumeration; however, 
can be decomposed into parts that are not invariant 
under permuting coordinates or spins Si separately. 
Let us first recap the interaction between a pair of mag- 
netic moments 
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rijSi 



Sj 'i{fij • Si){Vij ■ Sj 



(3) 



where 7 = geh/2mc (c is the velocity of light and the 
5-factor for dot carriers which even in few electron quan- 
tum dots is found to take basically bulk values). ^'^"^ Its 
Heisenberg-like first part is manifestly SU(2)-invariant 
in spin space and commutes with S^. This part nei- 
ther changes S nor Sz and just renormalizes the ener- 
gies slightly. It therefore can be ignored in view of the 
smallness of dipolar energies compared to the dot level 
separations. In Eq. ^ we abbreviate rij :— 



Ti — Tj and 
The second part of H]^^ can be decomposed 



-g-{rij ■ Si){rij ■ Sj) — 

ij 

where the three terms 

|2 



rriO) _ Iftjl 
rr(l) ^ C»i 



^S'+iS'-j -l- S'-iS'+j^ + C^ijSziSz 



Q—iji^Sj^iSzj + SziS^j) 



(4) 

(5) 
(6) 



Q+ij{S—iSzj + SziS-j 



rr(2) _ 



S+iS+j-, 
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are responsible to alter symmetries and spins after car- 
rying out summation over {i 7^ j). In Eqs. j?!) they 
change Sz by 0, ±1, and ±2, respectively. S± := SxUSy 
denote usual rising or lowering operators in spin space, 
g^ij := {xij ± ^yij)/r^j'^ is a complex (angular momen- 
tum generating) coordinate in the plane perpendicular to 
the axes of spin quantization, taken as the z-axes, and 

The spin changing part Eqs. ijHHZIl of can now 
further be decomposed according to partitions F of the 
symmetric group Sn, 



^ ^T=[N\ _^ ^r=[W-l,l] _^ ^r=[7V-2,2] 



(8) 



This latter representation is particularly useful to deduce 
non-zero transition matrix elements between quantum 
dot eigenstates of different total spins. No other par- 
titions occur since H^^ transforms as a product of two 
vector operators in position as well as in spin space, cf. 
Eq. (jSJl, i.e. as a tensor of rank two. In Eq. © i/'"=[^-i>i] 
changes the total spin S of dot by ±1 and i/'"=[^~2,2] 
±2, where the latter occurs only for A*" > 4 while the for- 
mer already for > 3. cannot achieve spin changes 
by more than ±2. For example, we can conclude already 
at this stage that the rate for direct transitions of an ex- 
cited S — Z quantum dot state into the (assumed) S = Q 
singlet ground state will be of the order 0{{H^Y) and 
therefore will be very small. All properly symmetrized 
operators H^=W-^^^\ and ffr=[A'-2,2] for TV = 3 and 
= 4 are listed in Appendix 1X1 

Note, that the property of H^^^^"^ to change Sz by 
0, ±1, ±2, respectively, is unrelated to their respective 
capability to change S. In the absence of further sym- 
metries of the quantum dot shape, all three operators 

contain both, i7r=[A'-i4] and i/r=[JV-2,2]_ 
cases of frozen electron motion in z-direction, as it applies 
to quantum dots (or arrays of quantum dots) fabricated 
on the basis of semiconducting hetero-structures^^ all of 
the above contributions involving vanish. Then H^^^^ 
simplifies and H'^^^ vanishes entirely, so that Sz can ei- 
ther remain unaltered (through i/'-^^) or change by ±2 
(through 



B. Two electrons 

Let us first focus on two electrons, i.e. N = 2. This is 
relevant, for example, for double dots containing one elec- 
tron on either side to realize the basic entity of coupled 
qubits,2*^ As already mentioned in the previous section, 
non-v4i symmetric partitions Eq. JSJ of occur only for 
N > 3. Further, does not contain the A2 partition 
for any N . Therefore, spin conversion transitions from a 
triplet excited state into the singlet ground stated will 
never be mediated by . In the related physics context 
of nuclear spin conversion of II2 molecules^ the stability 
of ortho-hydrogen (even over weeks) is traced back^ to 
parity symmetry of both, the molecule and the magnetic 
dipolar interactions between the two protons (of actually 
close proximity which enhances dipolar forces) to prevent 
the transition from the odd-parity ortho S — I state into 
the (by 80 Kelvin lower) even-parity para 5 = ground 
state. In this case, spins refer to the protons. In the con- 
text of quantum dots we can generalize this finding: Ir- 
respective of the shape of the quantum dot confining po- 
tential and of the functional form of the electron-electron 
interaction w(ri — r2) the dipolar interaction will not 
change (triplet or singlet) spin states as a result of per- 
mutational symmetry and quantum mechanical particle 
identity of iV = 2 electrons. This statement is not re- 
stricted to the lowest (golden rule) order 0{{H^)^) but 
even holds true to any order of . As one neat corollary 
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we conclude that two electrons in square shaped quantum 
dots (in the absence of other magnetic mechanisms) will 
stay in their respective spin states. This supports a cor- 
responding proposal for quantum computations based on 
superpositions of states where the two electrons occupy 
either of the two equivalent electrostatic energy minimum 
positions at diagonally opposite corners in a square4fi 

C. Strong interaction, pocket states 

In case of more than two electrons we focus on strong 
Coulomb forces at low carrier densities, i.e. at large val- 
ues of the electron gas parameter Tg ^ 1. Then, the 
kinetic energy is small and the electron system lowers 
its energy by Wigner localizing'^" the charge density near 
electrostatically favorable places. Precursors of Wigner 
crystallization have been found already at rg « 4 in 
two-dimensional quantum dots.^3_ A similar localization 
of charge density arises when the external confining po- 
tential separates the electrons, such as in the case of N 
quantum dots, each containing a single electron. In either 
case, at strong Coulomb interactions, eigenstates \4>n> 
of Hq Eq. (P) are well described by pocket statesrS42L2& 
which exploit the electron localization. They allow to 
estimate the spin dependent low energy spectrum to ex- 
ponential accuracy with increasing Tg, or with increasing 
dot separation. 

In the Wigner crystal state, electrons vibrate about 
electrostatic energy minimum positions. Linearizing the 
(Coulomb and external) forces yields the plasmon spec- 
trum of the confined A'^-electron system. Energy level 
separations ~ Wq -I- Ar~^ can be estimated from the 
dynamical matrix with a prefactor A depending on N 
and on the dot lay-out; is the confining frequency of 
the quantum dot(s). Due to the electron spin each plas- 
mon level is 2^-fold degenerate. Quantum corrections 
(partly) split this degeneracy into sub-levels, with all ex- 
hibiting well defined total spins S — 1 1/2 | : • ■ ■ j for 

{ odd } ^ [of + Zeeman de generacy, by spin 

rotation invariance], according to Sect. ^1 The ground 
state (in more than one spatial dimension) need not be of 
minimal spin Sq ~ or Sq ^ l/2]^^i^^i^^i^'^ A given spin S 
may appear more than once in such a spin split plasmon 
level; examples of spectra are discussed in Refs. 26..27,.28- 
The splitting arises due to permutational electron ex- 
changes by quantum mechanical tunneling through the 
electrostatic barrier consisting of the w-term plus the 
w-term in Eq. In the simplest case there are Nl 

different, but all energetically precisely equivalent, pos- 
sibilities to arrange the localized electrons; this defines 
I < p < Nl pocket states \p>. The width of each 
pocket state corresponds to plasmonic zero point oscil- 



lations and scales roughly as lu^^ ' in iVd-dimensional 
configuration space {d being the spatial dimensionality 
of the quantum dot, oftenS d = 2 but also d = 1 is 
realized, for example in rods of carbon nanotube^i). 
The energy scale A for spin splittings of plasmon lev- 
els through quantum mechanical electron exchanges is 
tuned by the magnitude of overlap integrals <p' \ Hq \ p> 
between two different arrangements p and p' . This latter 
quantity can be estimated semiclassicalljiSLSSi^i to read 
A ~ <p'\Ho\p> ~ ujpi exp(— y^) so that A/wpi < 1. 
Numerically obtained quantum dot spectra^SiSLi^ indeed 
nicely follow this behavior. For example, it exhibits 
the predicted^^ crossover into a spin polarized S" = 3/2 
ground state in a spherical two dimensional quantum dot 
containing N = 3 electrons at sufficiently low electron 
densityi^S'S^ As a result, all eigenstates 

|0n> = ^^c„p|p> (9) 

p 

belonging to the plasmon ground multiplet can ap- 
proximately be expressed through the set {|p>}. The 
(real) coefficients c„p, appearing in Eq. ensue from 
the irreducible representation F of the permutational 
group Sn, according to the wave functions symme- 
try which at the same time fixes the total spin S of 

I (/)„>; Nn = ^JY^pp' CnpCnp'<p'\p> cusures normahzation, 

<4>n\(l}n> = 1- 

D. Dipolar matrix elements 

Pocket states allow to conveniently estimate the ma- 
trix elements <0„ \H^\(j)n'> of the dipolar energy since, 
to leading order, electron positions may be taken as 
being well localized, i5-function like, for g^j and Qj in 
Eqs. [TJ or in the already symmetrized expressions in 
Appendix 1X1 This leads to a finite lattice spin problem. 
Having constructed symmetrized spin states, the matrix 
elements of for > 4 follow straightforwardly from 
Eq. ©. 

We demonstrate our approach for the particularly sym- 
metric cases of A^ = 3 and A = 4 electrons occupy- 
ing equilateral electrostatic equilibrium positions, as in 
a two-dimensional spherical quantum dot,^'^'^'^ in trian- 
gularly or square shaped quantum dots,"^^ or in equilat- 
eral triangular or square arrangements of single electron 
quantum dots. We assume frozen motion in z-direction, 
as in heterostructures, so that terms involving Qj or (f^ 
are irrelevant in Eqs. ((S}— |ZIl- Symmetrized, non-trivial 
spin states of minimal S'^-components are presented in 
Table D 
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iV 


5* 


Sz 


index 




■i 


3 
2 


1 

2 


A 

A 


:y5(.|TTi> + |TiT> + |iTT>J 
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1 
2 


1 
2 


Ea 


( 1 TTi> + e'^'"''^ 1 TiT> + e"'^''/^ 1 iTT>) 




3 


2 


2 


Eb 


71 ( 1 TTi> + e"'^'"/^ 1 TiT> + c'^''/^ UTT>) 




4 


2 





A 


^(|TTii> + UTn> + UiTT> + |TiiT> 4 


- |TiTi>+ UnT>) 


4 


1 





Tx 


5(lTTii> -iUTn> - UiTT> + i|niT>) 




4 


1 





Ty 


5(|TTU> + iUTn> - UiTT> - i|TUT>) 




4 


1 







;^(|TiTi> - UTiT>) 




4 








El 


5(|TTU> - UTTi>+ \Un>- |TiiT>) 




4 








E2 


;7|(|TTii> + UTn> + UiTT> + |TiJ.T> - 


-2|Tin>-2UTiT>) 



TABLE I: Symmetrized, non-trivial spin states of minimal ^z-component for TV = 3 and A'^ = 4. Spin states of larger \Sz\ are 
obtained easily. 



N ^3 



In this case X^i^j ofj ~ a.nd, in Appendix EI we 
replace Igy P by for = 3 where r is the mean inter- 
electron separation. Then, the only non-vanishing term 
takes the value 



^ 

3 ' 



rr(2) 



h.c. where "dij — 



is the azimuthal angle of Vij. 



Thus, necessarily changes Sz by ±2 and has non- 
vanishing matrix elements only between the A and the 
E states {5* = 3/2, & = 3/2} ^ {5 1/2, = -1/2} 
and {S = 3/2, Sz = -3/2} ^ {S = 1/2, S^ = +1/2} of 
Table n Their value emerges as: 



<0{A,S,=±3/2} 



In particular, this means that the not Zeeman aligned 
states of 5* = 3/2 with Sz — ±1/2 remain unaffected 
from dipolar decay. 



2. A = 4 

For a square arrangement of = 4 electrons two dis- 
tances occur: r along one edge and \/2r across the di- 
agonal. Inspecting aU the terms ijt'')!^^!] . . .ij(2)[2,2] j^j. 
A^ = 4 in Appendix ^ reveals that only contributions 



16 Sj5^J ^ 



non- vanishing. 



S. 



iS—j i S_^S_j 



h.c. remain 



term 
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Again, 5*^ has to change by ±2. The 
vanishes across the diagonal of the 



2 

0+ij - 



square, while g^^j — g^jj = along any edge. An eval- 
uation yields the non-zero matrix elements between sym- 
metrized spin states of Tabled i.e. 



<(t>{A,S^ = ±2}\H |0{T,,S.=O}> 
<0{A,S. = ±1} l'/'{T,,S,=Tl}> 



V2-32 

111 
64 r3 



<'p{A,S,=±2} 



{Ei,S,=0}> - ~lT7;-3 



If 

32 r3 



All other matrix elements vanish identically. In particu- 
lar, states of symmetries Tx, Ty, E2, and A with Sz = 
do not exhibit dipolar decay. 



E. Mixed spin states 

Non- vanishing elements of slightly mix eigenstates 
1 0„> of Ho- Here we are interested in admixtures to spin 
states S 



EE 

S'^S n 



> 



En,, — E„ 



(10) 

arising from other spins 5" ^ S. In Eq. (|10|l we have dis- 
regarded the very unlikely case of accidental degeneracy 
between eigenlevels of Hq (cf. (^) of different spins ^^"^ 

Eventually, this mixing will cause spin changing tran- 
sitions and thus spin relaxation. We disregard dipolar 
admixtures from other states | (f>n' > oi the same spin in 
Eq. pUfl as those occur much more efficiently by phonons, 
see in the subsequent Section. To this end, we take 
{|V'ng>} as exact eigenstates of Hq + H^. From now 
on we denote by rig the subset of n-values enumerating 
eigenstates of Hq that belong to the definite spin S. 



IV. SPIN RELAXATION RATES 

Typically, the electron-phonon interaction i/^i-ph es- 
tablishes thermal equilibrium between electron and lat- 
tice reservoirs on time-scales short compared to the times 
on which spin changing transitions occur. This is so be- 
cause the latter cannot be achieved directly by H'^^~p^ 
(cf. Sect. IIV'X|l . so that equilibrium will be established 
rapidly only among dot levels of given total spins. This 
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suggests to divide the total Hilbert space: 



of coupled elcctron-phonon states into orthogonal sub- 
spaces Ti-s, labeled according to the electron spin S. 
Transitions among subspaces occur only slowly by the 
action of while thermal equilibrium resides within 
each of the subspaces after much shorter times r°'^P^ 
at the lattice temperature (/cb/3)~^- Consider a certain 
electronic spin state S", as it may have been prepared, 
for example, using electronic transport techniques^*^^*^ 
Then, the rate 



(11) 



for its decay into a particular spin S" 7^ S" is given as 
the temporal increase of the spin S'-population {Ps{t))s', 
assuming an initial (i.e. after intra- Tig' equilibration has 
taken place) S" thermal equilibrium state, 

Ps'e-^"Ps'/TT{Ps^e-^"Ps'} . 

Here, Ps = Y.n^ I'/'n^X^rigl 1^'' projects onto Hs, 
IP^ denotes a unit operator on the phonon space. Tran- 
sition rates Rs^s' observe the detailed balance condi- 
tion, ensuring one vanishing (stationary) eigenvalue of 
the matrix Mss' = Rs'^S — Tg,^6ss' , which governs the 
rate dynamics. In the present context we are primarily 
interested in the total decay rate of the initial spin S' 
population; i.e., 



S^S' 



(12) 



In Eq. I|ll() the time evolution refers to the Hamiltonian 
H = i/o + i^° + -ffP^ + i?°'"P'' where the electron-phonon 
interaction i/°'~ph will be discussed next. This approach, 
in principle, accounts for rapid thermalizing spin conserv- 
ing multi-phonon transitions within subspaces Ti.s- 



A. Coupling to phonons 



wave lengths 2TTCs/ujp\ = 50 nm or 2ttCs/A = 500 nm, 
assumin g^^ o^p i ~ 3 meV and A ~ 0.3 meV, respectively, 
cf. Sect, nil CI for the definitions of the energies A and 
Wpi. At this wave lengths intra- valley scattering domi- 
nates. Its strength 



9l = 



El 



(14) 



mainly is regulated by the deformation potential con- 
stant E2 for longitudinal coupling which takes values of 
about^i^ 10 eV in Si. Further, g"^ depends on the mass 
density pM, the normalization volume V for the phonon 
modes, and on the sound velocity Cg. In Eq. H13|l the 
operator 



piq) - N 



(q) I (/-n^ >«/ 



(15) 



of the total electron density may excite the correlated 
electron system at non-zero q, though at conserved total 
spin S (and conserved z-component S^)- It can be de- 
composed into the basis {|(/!'n>}, where |0„> and \(j>n'> 
have same spin S, i.e., 

Pn^n'^,{q) = y^^^"^''/ dra . . .drAr<r,r2, . . . ,rAr|0„^> 

\r,r2, ■ ■ ■ ,rN> . (16) 



X<(7 



At small q 

Pnn'{q) = <5, 



E \q\ these coefficients are expanded, 
-|- a{qHY where, to lowest non- vanishing 
order, v — 1 unless the electron charge density distribu- 
tion of the quantum dot or of the ensemble of quantum 
dots is parity symmetric, in which case v = 2. The quan- 
tity i either equals the typical distance between electrons 
if n and n' belong to the same plasmon multiplet, or 
t ~ (mwpi)"^/^ for n and n' from different plasmon mul- 
tiplcts. The magnitude of a can be estimated by inserting 
Eq. © into Eq. (|l(j|l and using, for convenience, the den- 
sity distribution Pnn'{f) iu real space. This reveals that 
a is proportional to the maximum overlap between un- 
equal pocket states, i.e. max <p I p'>, a quantity which, 

in turn, is proportional^ to A/wpi. 



Electron-phonon coupling in semiconductors has been 
studied intensively in the 1950-ies and 1960-ies. For ho- 
mopolar semiconductors, such as Si or Ge, deformation 
potential coupling^ has been established. It can be ex- 
pressed as^'' 

H'^^-P^ = J29iPiQ)(b, + bt,) (13) 
q 

where we have suppressed the phonon branch index. 
Considerably below room temperature, pertinent to pos- 
sible quantum computing, optical phonons don't con- 
tribute so that in Eq. (|13|l is meant to create a lon- 
gitudinal acoustical phonon of momentum q. For excita- 
tions of the electronic system most relevant are phonon 



B. Transition rates 

We are now in the position to calculate the phonon 
mediated transition rate Eq. (|lll) as a result of spin mix- 
ing, see Eq. H1U|) . Assuming a not too strong electron- 
phonon coupling, use of standard time dependent per- 
turbation theory with respect to TJ^i-ph^ explicated in 
Appendix IbI yields to leading order the rate Rs^s', cf. 
Eq. Hll|) . reading 

Rs^s' = f; E e"'^^".'J„^,„,(|£;„^,-i?„J) 

X K|£^n^, - ^«J) + e(i?„^, - i?„ J] , (17) 
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where we have defined the (temperature independent) 
coupled density of phonon states for transitions between 
spins S and S", respectively, as 



Jn^,n^{uj) = Xl^'q ^(^ " l^l^s) 



En,„ — En' 



(18) 



In Eq. ifTTjl Zs = J2n^ ^ ^^"s denotes the partition 
function inside the subspace Tis, n{uj) = {e^^ — 1)^^ 
the Bose function, and 8(x) the Heavyside step func- 
tion. Straightforwardly, higher order terms regarding 
^ci-ph g^jgQ considered for Rs^s'i although the 
corresponding explicit expressions are rather lengthy. In 
Eq. pUI) we have assumed that phonon states and en- 
ergy eigenvalues remain unaffected by the weak dipolar 
mixing. 

At low temperatures, T <C A/fce compared to the typ- 
ical distance A between dot levels of same or of differ- 
ent total spins inside the lowest plasmon multiplet, the 
Bose factor n(|A|) ^ 1 is small and only the ground 
level rigg will be occupied within each subspace Tis- In 
this temperature regime thermalization into the global 
ground state n — Og^ of spin Sq will take place ex- 
clusively through the direct process by emission of a 
resonant phonon of energy A so that the coupled den- 
sity of states J„ (^) controls the relaxation rate 
Rso^s- Still, a summation over excited levels n'g^ > 
and riQg > n^g appears in Eq. H18|l . as the lowest 
terms^^ = s-nd n'^g = n^g cancel exactly. In Sec- 
tion IIV a] it has been estimated that non-diagonal co- 
efficients pnn'iq) (A/tt'pi)(rA/cs)'' for low energy and 
long wave lengths transitions; here r denotes the distance 
between electrons and = 1 or = 2 in the absence or 
presence of parity symmetry. For the electron-phonon 
coupling Eq. (|14|l this results in a spin transition rate at 
zero temperature through the direct process, reading 



RSo^ 



— 27r J„ 



,Os„(A) 



(19) 



with 



In Si the rate Eq. (|19(l appears to be very small at 
zero temperature, ~ 10"^ s"^ for three electrons at den- 
sities corresponding to rg — 1, and considering a quan- 
tum dot'^^ of LOp\ ~ 3 meV and A = 0.3 meV. However, 
this number strongly varies with parameters as seen in 
Eq. (|19|l . Parity symmetric quantum dots (where v — 2) 
would suppress this decay rate even further at small tran- 
sition energies due to J„o(A) ~ A^ in this case. These 
values are, of course, considerable smaller than the decay 
rates estimated from spin orbit effectsjiSiS if present. 

They are also smaller than the rates estimated from the 
hyperfine interaction with nuclei of non-zero spin^^ti^ii^ 
Particularly in Ref. ^3 a hybrid mechanism is consid- 
ered which is closely related to the one presented here 
in combining the electron-phonon coupling with a spin- 
mixing interaction. Transitions between total spin 5 = 1 
and S* = of a two-electron quantum dot are investi- 
gated. The low temperature rate has been estimatedi^ 
to ~ 10~^ s~^ for similar quantum dot parameters as 
above, assuming a two-dimensional dot fabricated on the 
basis of heterostructures. This rate is proportional to the 
number N-^ of non- vanishing nuclear spins covered by the 
electron wave function. In GaAs almost every nucleus 
has spin / = 3/2. It is instructive to determine from this 
result a critical concentration of ^^Si nuclei in sili- 
con, the only ones of non- vanishing spin / ~ 1/2, beyond 
which the here described dipolar mechanism should pre- 
vail over the spin decay via nuclear spins. For a quantum 
dot of the same excitation energy, the electron wave func- 
tion in natural silicon covers only about A^f ' « 10^ of the 
29Si nuclei whilei^ in GaAs Af^^^" w 10^ Two further 
important differences between Si and GaAs have to be 
taken into account. Firstly, the type of electron-phonon 
coupling which is piezo-elastic in GaAs while we have de- 
formation potential coupling in Si. Accidentally, for the 
here considered quantum dot parameters (and assuming 
again heterostructures and now laterally parabolic con- 
fining potential) J„„'(u;) in Si is only by 0.6 smaller than 
in GaAs. Secondly, the nuclear spin — 1/2 of ^^Si as 
compared to /^aAs _ 3^2 in GaAs which reduces the cou- 
pling by /Si(/Si -I- i)//GaAs(-jGaAs + 1) = 1/5. This yields 
w 2 X lO"'^ nm-3 (note that Cf w 2.5 x 10"^ ^m-^ 
has been reported^" experimentally). This value is less 
stringent than the isotopic purification required for the 
quantum computer^^ based on the nuclear spins of ^^P 
donors, where Cn should be smaller than 10"'' nm^"^ 
in Si with N being the number of qubits. 



Eh' 



C. Temperature Dependence 



unless this transition is not suppressed entirely for cases 
discussed in Sect. ImDl In Eq. lO we have assumed 
for simplicity that level separations En' — Eq„ ^ A 

Sq 

and — En ~ A both are of the order— A. Also, 

"s OS ' 

we have inserted the areal density ria = ''"^ of electrons, 
focusing on the measured quantity in two-dimensional 
samples. 



Through the marked increase of J„o(t^) 



as a 



function of transition energy w, relaxation can take ad- 
vantage from spin conserving thermal excursions to plas- 
monic excited levels and accomplish the spin transition 
at an elevated energy. In NMR-theory this possibility is 
called the 'Orbach process— and shows up in a steeply 
increasing relaxation rate with temperature. Our formu- 
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lation, Eq. H17() . of the transition rate explicitly incorpo- 
rates such thermal excursions. They turn out to influ- 
ence considerably the temperature dependence of Rso^s 
which only at low temperatures follows the Bose-behavior 
^ J„o(A)n(A) of direct transitions. Already at tempera- 
tures not much exceeding A the severely stronger increase 
~ Jno(^pi) exp(— Wpi/fcBT), following from Eq. (|17(l . can 
easily enhance the transition rate by three orders of mag- 
nitudes, depending on system parameters. A similarly 
pronounced increase of (nuclear) spin relaxation rates 
has been discussed in detail^ in the context of quantum 
rotating molecules: substantial increases in spin chang- 
ing transition rates by more than six orders of magni- 
tudes are depicted with Figure 3b of Ref. Again, the 
stronger increasing density of coupled phonon states in 
quantum dots of parity symmetry should lead to even 
more pronounced temperature sensitivity. 

V. RESUME 

We have investigated dipolar interactions between the 
magnetic moments of electrons confined to one or to sev- 
eral quantum dots and studied the rate of inelastic to- 
tal spin changing transitions. As compared to the cou- 
pling to nuclear spins lii^iiSiiSiS and to spin-orbit in- 
duced decay,^'^'"'^^'^^'^- dipolar spin decay turns out as 
much weaker. However, either of the earlier studied 
mechanisms can, at least in principle, be eliminated by a 
suitable sample design. It is therefore possible that the 
dipolar interactions between electronic spin moments, to- 
gether with the coupling to lattice modes, as discussed 
in the present work, will ultimately limit long time quan- 
tum computations, even when devices become optimally 
designed. Experimentally, the dipolar mechanism should 
show up most directly by observing at low temperatures 
the dependence on the electron density, cf. Eq. l(T^ . 

Upon generalizing previous approaches we incorporate 
here electron-electron interactionsiiii^ in the quantum 
dot(s) which, additionally, are important for magnetic 
featuresiS2i2iS4-25-3i Pqj- example, ground state spins may 
exceed the values 5o = or 5*0 = 1/2 expected for non- 
interacting electrons. We focus on the limit of strong 
interactions, where electronic many-body wave functions 
can be described as 'pocket states i^StSLS^ and where the 
spectrum exhibits spin-split plasmon multiplets. The 
dipolar interaction is decomposed according to the sym- 
metric group and non- vanishing matrix elements are de- 
termined in their dependence on spatial parts of the col- 
lective electron wave functions. As we have shown, par- 
ticular spatial symmetries of the quantum dot(s) can re- 
duce the number of non-zero elements. 

One important result is the stability of = 2 elec- 
tron spins. Irrespective of the dimensionality, the shape 
of the quantum dot(s), or of the electron-electron in- 
teraction strength, the decay of triplet states into the 
singlet ground state is always suppressed. Any dipolar 



decay channel will require participation of further elec- 
trons. This is important, for example, for two coupled 
qubits, the basic (gate)-element for quantum information 
processing. 

Also at larger electron numbers, non-ground state 
spins can be stable with respect to dipolar interactions. 
We have discussed the case of = 3 electrons on an 
equilateral triangle. Here, Sz = ±1/2 states of the 
5 = 3/2 sub- manifold prove robust against decay into 
any 5* = 1/2 state. Only spin polarized Sz = ±3/2 states 
decay into the 5=1/2 sub-manifold. Further, A'' = 4 
electrons on a square shaped quantum dot exhibit robust 
5 = 1 states of Tx and Ty symmetry, and 5 = states 
of E2 symmetry. 

Owing to the smallness of magnetic and in particular 
dipolar energies, compared to dot level separations, the 
energy accompanied with an actual spin transition has 
to be provided by the reservoir of lattice vibrations. As 
in previous approachesi'S'SiiliiH we have considered the 
coupling to acoustic phonons. Parity symmetric dots are 
weaker coupled to phonons, which further suppresses spin 
decay in this case. Additionally, we have accounted for 
rapid thermal excursions of the system within electron- 
phonon subspaces of given (many electron) dot spins. 
This enables one to deduce the dependence of spin relax- 
ation over a wider range of temperatures as compared to 
the resonant direct process. As a result we found a very 
striking increase of the spin decay rate. This rate grows 
with temperature considerably steeper than the naively 
expected proportionality to the Bose function describing 
direct processes: It occurs already at temperatures that 
barely exceed the energy difference between the lowest 
levels of different spins, but is still considerably smaller 
than the energy for plasmon excitations. Although we 
find amazingly stable spin configurations at low temper- 
atures this marked temperature sensitivity restricts the 
operation temperatures of quantum computing dots (un- 
less quantum computation can be confined to the sta- 
ble spin configurations) to values that are not exceeding 
much the lowest level separations. 

Because the same phonon energy reservoir is consid- 
ered in previous workSiSiSiiliiiii^ for spin decay, a similar 
scenario regarding the dot symmetries and the spin con- 
serving phonon induced excursions should apply also to 
magnetic mechanisms of the spin-orbit or of the hyperfine 
type. We expect therefore a similarly striking tempera- 
ture sensitivity as obtained here for these mechanisms as 
well. 
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APPENDIX A: SYMMETRIZED OPERATORS INDUCING SPIN CHANGING TRANSITIONS FOR 

N ^3 AND iV = 4 ELECTRONS 
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The following explicit form of dipolar operators Eq. © contain non- vanishing elements for A'' = 3: 



7J(0)[2,1] 



ij(i)[24] 



ij(2)[2,l] 



(2) _ _1_ 

y 24 



EC 



^(5'+iS'_j + S-iS+j) 



^^(S'+iS'zj + SziS^j) 



E Ci t'+^J 



E 



J2 



E ^'-y 



E 
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and for iV = 4: 

_H-(0)[3,1] 



ij(2)[3,l] 
ij(0)[2,2] 



+ [Cij9+ij ~ Cijfi'+ij] {S-iSzj + SziS-j) — {S_iS^j + S^iS^j) 



5* iS j S r^S j 



2l 



_H-(1)[2,2] 



12' 



EC 



+iSzj + SziS+j) + {S^iSzj + SziS^j) 



+ [Cii£'+ij + Cyfi'+ij] {S-iSzj + SziS-^j) + {S_iS^j + S2iS_j) 



24' 



^^(S'+i'S'zj + SziS+j) 



^ ' ^ ^ {S—jSzj + SziS-j 



_H-(2)[2,2] ^ 



1 

16 



E^ + 1 I '5'-«5'-j + S_jS_j 



48' 



E ^1^1 



E ^-■'^-j 



In the above expression, (i, j) take the two values out of 1, ... ,4 that are both different from (i, j). 



APPENDIX B: DERIVATION OF EQ. ([Tft 



In a perturbative expansion w.r.t. H'^^ p'' of either of the two time evolution operators appearing in Eq. Hll|l we 
write 



lP''-i / dt' e'^o^'lf^'-P'^e-'^o*' 
Jo 



(Bl) 



t ft' 

At' \ dt" e'^''*'iJ°'"P'^e"'^''(*'"*"^i7'''~P''e"'^»*" 
Jo 
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Here, Ho = Hq + + H^^ with Hq defined in Eq. (QJ and the eigenstates of Hq + taken according to Eq. ifTUIl . 
To second (i.e. lowest non-vanishing) order in iJ^'^P^' only two of the second terms in the square bracket of HB1() 
contribute to Eq. ((TT|l . yielding with (fT!^ 



1 V- 2 -/^^...'^ / dt' I dt"i I e" 



x(-i) <0„ le'(^°+«°)*"p(q)e-'(^°+^°)*" |0„,,> ke-l-^K*'-*") + (n, + l)e--=l'^l(*'-*") 



Here, the Bose factors fiq = (e'^'^='''l — 1)~^ result after thermal averaging over phonon modes. Inserting now eigenstates 
\il)n> of {Hq + H^) for the |0„>, according to ((Tn|l . and carrying out the long time hmit t oo, yields for the decay 
rate 



Rs^S' 



2tt 



i <^n„, I - 



2-^ E — F, 

ngii S s 



X |V'ns> 



<0ng„ I g° I (j)ns > 

^ E — E 

n ,, "s" 



-i <i'ns I - 



B — B 



\ "S" 



X [nq5(Ens, - Ens + Cs|g|) + (f^q + l)<5(£^ns, - £^ns " Cs|q|)] , 

which is readily brought into the form p7|l with ^18B . after employing p5|) and using the orthogonality of spin states. 
Further calculational details can be found in Ref. js^- | 
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